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We report results on a Dalitz analysis of three-body charmless B^ —> decay including 

searches for direct CP violation. We report the first observation of the decay B^ —> /2(1270)it'^ 
with a statistical significance above 6 ct. We also observe 3.9 ct evidence for large direct CP violation 
in the B^ —> p{770)° channel. This is the hrst evidence for CP violation in a charged meson 
decay. The results are obtained with a data sample that contains 386 million BB pairs collected at 
the T(45') resonance with the Belle detector at the KEKB asymmetric-energy e"'"e“ collider. 

PACS numbers: 13.25.Hw, 11.30.Er, 14.40.Nd 


Decays of B mesons to three-body charmless hadronic 
final states provide new possibilities for CP violation 
searches. In decays to two-body final states {B —> Ktt, 
TTTT, etc.) direct CP violation can be observed as a dif¬ 
ference in B and B decay rates. In decays to three-body 
final states that are often dominated by quasi-two-body 
channels, direct CP violation can also manifest itself as a 
difference in relative phase between two quasi-two-body 
amplitudes that can be measured via amplitude (Dalitz) 
analysis. So far direct CP violation has been observed 
only in decays of neutral K mesons P| and recently in 
neutral B meson decays 0- However, large direct CP 
violation is expected in charged B decays to some quasi- 
two-body charmless hadronic modes 0 . 

The search for direct CP violation in the three-body 
charmless B^^K^tt^tt^ decay described in this Let¬ 
ter is performed by applying a Dalitz analysis tech¬ 
nique y to a data sample containing 386 million BB 
pairs, collected with the Belle detector operating at 
the KEKB asymmetric-energy e“'"e“ collider 0 with a 
center-of-mass (c.m.) energy at the T(4S') resonance. 
These results supersede the results reported in Ref. 0. 

Charged tracks are required to have momenta trans¬ 
verse to the beam greater than 0.1 GeV/c and to be con¬ 
sistent with originating from the interaction region. For 
charged kaon identihcation we impose a requirement on 


a particle identification variable which has 86% efficiency 
and a 7% fake rate from misidentihed pions as measured 
from data. Charged tracks that are positively identihed 
as electrons or protons are excluded. B candidates are 
identified using two kinematic variables: the energy dif¬ 
ference AE = Vc^lPip + Cml) - and the 

beam constrained mass Mbc = ^ SiPiP) 

where the summation is over all particles from a B can¬ 
didate; Pi and rrii are their c.m. three-momenta and 
masses, respectively; I-beam i® beam energy in the 
c.m. frame. The signal Mbc resolution is mainly given 
by the beam energy spread, and amounts to 2.9 MeV/c^. 
The signal AE shape is fitted by a sum of two Gaussian 
functions with a common mean. In hts to the experimen¬ 
tal data, we hx the width (35 MeV) and the relative frac¬ 
tion (0.16) of the second Gaussian function to the values 
obtained from Monte Carlo (MG) simulation. The com¬ 
mon mean of the two Gaussian functions and the width 
of the main Gaussian are allowed to float. 

The dominant background is due to e^e~^qq {q = 
u, d, s and c quarks) continuum events. We reject about 
98% of this background while retaining 36% of the sig¬ 
nal using variables that characterize the event topology. 
For more details see Ref. 0 and references therein. From 
MG studies we find that the dominant backgrounds orig¬ 
inating from other B decays that peak in the signal re- 
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FIG. 1: AE distributions for (a) B~ ^K~n~-K^ and (b) 
events with |Mbc—Ms[<7.5 MeV/c^. Points 
with error bars are data; the smooth curve is the ht result; 
the hatched areas are various background components. 


gion are due to ] 7 r+ and due to 

in which muons are misidenti- 
fied as pions. We veto these backgrounds by apply¬ 
ing requirements on the invariant masses of the appro¬ 
priate two-particle combinations 0. The most signif¬ 
icant backgrounds from charmless B decays originate 
from B'^—> ri'['yTT'^'K~]K'^^ B~^^ tt~ where one of 
the two same-charge pions is misidentified as a kaon, and 
from B^^K'^'k~ processes. 


The AE distributions for B^ TT^TT^ candidates 

that pass all the selection requirements are shown in 
Fig.ca In the fit to the AE distribution we fix the shape 
of the BB background component from MC and let the 
normalization float. The shape of the qq background is 
parametrized by a linear function with the slope and nor¬ 
malization as free parameters of the fit. From these fits 
we find 2248 ±79 (2038 ±76) B~ {B'^) signal events; the 
width of the main signal Gaussian is 15.3 ± 0.5 MeV. 

For the amplitude analysis we select events from the 
B signal region defined as an ellipse around the Mbc and 

AF; mean values: ±[ 4 cniv]"<l- The to¬ 

tal number of events in the signal region is 7757; the rel¬ 
ative fraction of signal events is 0.512 ±0.012. The distri¬ 
bution of background events is determined from analysis 
of events in the M^c—AE sideband region. 

The analysis is performed by means of an unbinned 
maximum likelihood fit. The distribution of background 
events is parametrized by an enmirical function with 11 
parameters 0 • As found in Ref. |3, the three-body > 
TT amplitude is well-described by a coherent sum 
of R:*(892)07r+, F:^(1430)°7r+, p{770)°K+, /o(980)F:+, 
/x(1300)iF+ and XcoK'^ quasi-two-body channels and a 
non-resonant amplitude. The /x(1300)iF+ channel was 
introduced in order to describe an excess of signal events 
at M( 7 r+ 7 r“) ~ 1.3 GeV/c^ (see Fig.l^b)). The best fit 
is achieved assuming /jf(1300) is a scalar state; the mass 
and width determined from the fit (see below) are con¬ 
sistent with those for /o(1370) na. Each quasi-two-body 
amplitude includes a Breit-Wigner function, a B decay 
form-factor parametrized in a single-pole approximation, 


a Blatt-Weisskopf factor for the intermediate reso¬ 
nance decay, and a function that describes angular corre¬ 
lations between final state particles. This is multiplied by 
a factor of that describes the relative magnitude and 
phase of the contribution. The non-resonant amplitude is 
parametrized by an empirical function = 

where a, a”’’ and Sf'^ are 
fit parameters, S 13 = and S 23 = M^( 7 r+ 7 r“). 

In this analysis we modify the model by changing the 
parameterization of the /o(980) lineshape from a Breit- 
Wigner function to a Flatte parameterization and by 
adding two more channels: tu(782)iF+ and /2(1270)R''^. 
For CP violation studies the amplitude for each quasi- 
two-body channel is modified from ae*^ to ae*‘^(l ± 6 e®‘^), 
where the plus (minus) sign corresponds to the R"*" {B~) 
decay. With such a parameterization the charge asymme¬ 
try, Acp, for a particular quasi-two-body B^f channel 
is given by 


N — iV+ 2bcosip 
= N- +N+ ^ T+W' 


( 1 ) 


It is worth noting that in this parametrization we assume 
zero relative phase between B~ and B'^ amplitudes. 

First we fit the data with bi = 0 (no CP violation) and 
determine the parameters of the /x(1300) (M = 1.449 ± 
0.013(stat.) GeV/c^ F = 0.126 ± 0.025(stat.) GeV/c^), 
/o(980) (M = 0.950 ± 0.009(stat.) GeV/c^ and cou¬ 
pling constants = 0.23 ± 0.05(stat.) and pkk = 
0.73±0.30(stat.)), and the parameter of the non-resonant 
amplitude a = 0.195 ± 0.018(stat.). We then fix these 
six parameters and repeat the fit to data with b and 
(p floating for all terms except B^^ u!{782)K^ and the 
non-resonant amplitudes. Possible effects of these as¬ 
sumptions were studied, and are included in the final 
results as a part of the model uncertainty. Projections 
of the fit are shown in Fig.[3 and the results are sum¬ 
marized in Tabled We find that the statistical signif¬ 
icance of the B^^f2{1270)K^ signal exceeds 6 cr; this 
is the first observation of this decay mode. The signifi¬ 
cance of the B^^ lo{782)K^ signal is 2.1tT. The statis¬ 
tical significance of these signals (and the asymmetries 
quoted in TableP) is calculated as a/— 21 n(£o/>Cmax), 
where /Imax and £0 denote the maximum likelihood with 
the nominal fit and with the corresponding amplitude 
(or asymmetry) fixed at zero, respectively. Note that 
the significance of the asymmetry is sensitive not only 
to bi cos ifi , but also to bi sin ipi , whereas Acp is sensi¬ 
tive only to the former (Eq.^l. Therefore, the signifi¬ 
cance of the asymmetry should be interpreted as hav¬ 
ing two degrees of freedom. The only channel where 
the significance of the asymmetry exceeds the 3a level is 
B^^p{770)°K^. Figures|2Ka,b) show the M( 7 r+ 7 r“) dis¬ 
tributions for the p(770)—/o(980) mass region separately 
for B~ and B'^ events. However, the interference term 
between B^p{770)K vector and B^fo{980)K scalar 
amplitudes cancels out when making the pro- 
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TABLE I: Results of the best fit to events in the B signal region. The first quoted error is statistical and the second 

is the model dependent uncertainty. The quoted CP asymmetry significance is statistical only. 


Channel 

Fraction (%) 

5(°) 

b 

<9 0 

Asymmetry significance (a) 

A'*(892)7r± 

A'o*(1430)7r± 

p(770)°A'± 

13.0 ± 0.8/^'J 
65.5±1.5/|^ 
7.85 ± 0 . 93 / 0:59 

0 (fixed) 

55 ± 4+1 
-21 ± 14/)g 

U.U/0_Q Q3]^_Q QQ3 

A A£^Q +0.032+0.010 

U.UUy_Q Q3Q_Q QQ3 

^ <^0+0.12+0.07 
U.zo_o.o9-0.09 

-18 ± 44/53 
-123 ± 16/^ 
-125±32/g^ 

2.6 

2.7 

3.9 

tj(782)A'± 

/o(980)A± 

/2(1270)A'± 

0.15±0.12/!(;(;^ 
17.7± 1.6/3;3 
1.52 ± 0.35/!(-3? 

100 ± 31 / 2 ? 
67± ll/CJ 

140 ± ll/J,® 

0 (fixed) 
0.30±0.19/(;;?g 

^ Qiy+0.19 + 0.11 
-0.16-0.04 

-82 ± 8/2 
-24 ±29/1^ 

1.6 

2.7 

/x(1300)A± 

4.14±0.81/“:^J 

-141 ± 10/9 

0.12±0.17/(J:S^ 

-77 ±56/13 

1.0 

Non-Res. 

34.0±2.2/i;8 

ST = -n±5tl 
= 185 ± 20/i9 

0 (fixed) 

— 

— 


TT 2 TM 2 ” 

-118 ±24/^' 

0.15 ± 0 . 35 / 0 , 0 ? 

-77±94/;f 

0.7 


jection for the entire range of the helicity angle of the 
7r“'"7r“ system (the angle between the kaon and the pion 
of the opposite charge in the tt+tt” rest frame.) Thus 
only the difference in relative fractions can be observed 
from comparison of Figs. Eta) andOtb). The effect is more 
apparent in spectra for the two helicity an¬ 

gle regions cos9'Jf<0 and cos9'^>0 shown in Figs.^tc- 
f). Here the difference in the interference terms for the 
B~ and 5+ decay amplitudes (due to different relative 



FIG. 2: Results of the fit to events in the B sig¬ 
nal region: (a) M{K^■k~) distribution with > 

1.5 GeV/c^; (b) distribution with > 

1.5 GeV/c^. Points with error bars are data, the open 
histogram is the fit result and hatched histogram is the 
background component. Inset in (a) shows the A* (892) — 
A^KldSO) mass region in 20 MeV/c^ bins. Inset in (b) shows 
the XcO mass region in 10 MeV jc? bins. 


phases between the B^p{770)K and B^fo{980)K am¬ 
plitudes) can be distinguished as a difference in the shape 
of the M(7r+7r“) spectra for B~ and decays. Re¬ 
sults of the branching fraction and Acp measurements 
are summarized in TablelTTl where for intermediate res¬ 
onance fractions we use world average values 0. The 
reconstruction efficiency is 22.4±0.2%, determined from 
signal MC simulation in which events are generated ac¬ 
cording to the matrix elements obtained from the best fit 
to data. 

To assess how well any given fit represents the data, 
the Dalitz plot is subdivided into non-equal bins requir¬ 
ing that the number of events in each bin exceeds 25. 
A pseudo-y^ variable for the multinomial distribution is 
then calculated as = — 2 > where Ui 

is the number of events observed in the f-th bin, and Pi is 
the number of predicted events from the fit. More details 
are given in Ref. 0. The x^/-^bins value of the fit to sig¬ 
nal events is 182.5/141 (32 fit parameters) and 127.6/120 
for the fit to background events. 

The following sources of systematic error are found 
to be dominant in the determination of branching frac¬ 
tions: charged track reconstruction (3% in total); parti¬ 
cle identification efficiency (4.5% in total); requirements 
on event shape variables (2.5%); signal yield determina¬ 
tion from the AE fit (3.9%); model dependence (1%); 
number of produced BB pairs (1%). For the quasi-two- 
body channels additional sources are the uncertainty in 
parametrization of the background density function and 
the uncertainty in secondary branching fractions (2% for 
/2(1270), 11% for A:*(1430) and 10.8% for XcO [III)- Note 
that in the asymmetry calculation most of these system¬ 
atic uncertainties cancel out. A few remaining sources are 
uncertainty due to a possible asymmetry in background 
from charmless B decays (0.6%); the possible bias due 
to intrinsic detector asymmetry (1.6%); B^ signal yields 
determination (1.1%). 

To estimate model uncertainty in the branching frac¬ 
tions and Acp for individual quasi-two-body channels. 
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TABLE II: Summary of branching fraction results. The first quoted error is statistical, the second is systematic and the third 
is the model uncertainty. Note that XcoK'^ contribution is not included in the three-body charmless branching fraction. 


Mode 

B(B± ^ ^ A±7r±7r=F) x 10® 

S(B± ^ Rh^) X 10® 

Acp (%) 

Charmless 

A*(892)[A±7r=F]7r± 

Ao*(1430)[A±7rT]7r± 

p(770)°[7r+7r“]A^ 

/o(980)[7r+7r-]A± 

/2(1270)[7r+7r-]A± 

Non-resonant 

48.8 ± 1.1 ± 3.6 

6.45 ± 0.43 ± 0.48tg;3t 

32.0 ± 1.0 ± 2.4ti;^ 

3.89 ± 0.47 ±0.29lK;^g 

8.78 ± 0.82 ±0.65t?;6i 

0.75 ± 0.17 ±0.06i;^;(8 

9.67 ± 0.64 ± 0.72to f2 
51.6 ±1.7 ±6.813 ? 
3.89 ± 0.47 ± 0.29l!);?g 

1.33 ± 0.30 ±0.111!);3^ 
16.9 ± 1.3 ± 1.3lo g 

±4.9 ±2.6 ±2.0 
-14.9 ± 6.4 ±2.0lK j 
±7.6±3.8±2.0lg;^ 
±30 ± 11 ± 2.01“ 
-7.7±6.5±2.0llg 
-59±22±2.0l? 

Xco [7r+7r“]A=^ 

0.56 ± 0.06 ± 0.041“;^^ 

112 ± 12 ± 18lf 

-6.5 ± 20 ± 2.01?;^ 


we vary the nominal model and repeat the fit to data. 
The following variations are performed separately: we 
add one additional channel which is either Ar*(1410)°7r+, 
Ar*(1680)'’7r+, or Ar|(1430)°7r+; remove w(782)A'+ or 
/2(1270)itr+ channel from the nominal model; fit the 
data assuming /jf (1300) is a vector (p(1450)) or exclud¬ 
ing this contribution; and use several alternative param- 
eterizations of the non-resonant amplitude ll 2l | . To 
cross check the asymmetry observed in B^^p{770)^K^, 
we make an independent fit to B~ and B^ subsam¬ 
ples. We also confirm the significance of the asymmetry 
observed in B^^p(770)^K^ channel with MC pseudo¬ 
experiments where events are distributed according to 
the matrix element determined from the fit to data. All 
the cross-checks give consistent results. Finally note that 



FIG. 3: 7r^7r“ mass spectra for B~ (left column) and B'^ 
(right column) events for different helicity regions: (a,b) no 
helicity cuts; (c,d) cosOTT^ < 0; (e,f) cos 6^ > 0; Points with 
error bars are data, the open histogram is the ht result and 
the hatched histogram is the background component. 


the second solution with a much smaller fraction of the 
Ar*(1430)7r signal as found in Ref. 0 is confirmed in this 
analysis. However, comparisons with results on elastic 
K-t: scattering 0 and with some theoretical considera¬ 
tions 0 favor the solution with a large iF*(1430)7r frac¬ 
tion. We find that values of CP parameters b and (p are 
almost solution-independent; variation in their values is 
considered as a part of the model uncertainty. 

In conclusion, we have performed an amplitude anal¬ 
ysis of the three-body charmless B^ decay. 

The branching fractions for a number of quasi-two-body 
channels have been measured; we report the first obser¬ 
vation of H+^/2(1270)iF+, a tensor-pseudoscalar decay. 
We also perform a search for direct CP violation in quasi- 
two-body intermediate states and find evidence for large 
direct CP violation in the decay B^^p{77Cf'. This 
is consistent with recent results from BaBar [l^ and with 
some theoretical predictions . The statistical signifi¬ 
cance of the asymmetry is 3.9tT and varies from 3.7cr to 
4(7 depending on the model used to fit the data. This 
is the first evidence for CP violation in the decay of a 
charged meson. Statistical significance of the Acp for 
the B^^p{77C)^K^ channel is 3.0(7. 
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